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Abstract 

For ordinary (independent) percolation on a large class of lattices 
it is well-known that below the critical percolation parameter pc, the 
cluster size distribution has exponential decay, and that power-law 
behaviour of this distribution can only occur at pc- This behaviour is 
often called 'sharpness of the percolation transition'. 

For theoretical reasons as well as motivated by applied research, 
there is an increasing interest in percolation models with (weak) de- 
pendencies. For instance, biologists and agricultural researchers have 
used (stationary distributions of) certain two-dimensional contact-like 
processes to model vegetation patterns in an arid landscape (see [19] ) . 
In that context, occupied clusters are interpreted as patches of vege- 
tation. For some of these models it is reported in [19] that computer 
simulations indicate power-law behaviour in some interval of positive 
length of a model parameter. This would mean that in these models 
the percolation transition is not sharp. 

This motivated us to investigate similar questions for the ordinary 
('basic') 2D contact process with parameter the infection rate A. We 
show, using techniques from the papers [8] and [11] by Bollobas and 
Riordan, that for the upper invariant measure v'x of this process the 
percolation transition is sharp: If A is such that (i^A-a.s.) there are no 
infinite clusters, then for all parameter values below A the cluster-size 
distribution has exponential decay. 
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1 Introduction and statement of the main re- 
sult 

The Contact Process was introduced as a stochastic model for the spread of 
an infection in a population with a geometric structure, usually represented 
by the (i-dimensional cubic lattice. Each vertex x of this lattice represents 
an individual, whose state, infected (1) or healthy (0), at time t is denoted 
by (Tx{t)- The dynamics in this model is as follows: A vertex in state 
goes to state 1 ('becomes infected') at a rate equal to A times the number 
of neighbours of that vertex that are in state 1. A vertex in state 1 goes to 
state ('recovers') at rate 1. Here A is the parameter of the model, called 
the infection rate. In this paper we restrict to the case d = 2. Depending 
on the applications one has in mind, the terms 'infected' and 'healthy' are 
sometimes replaced by 'occupied' and 'vacant' respectively. In the remainder 
of this paper we will use this latter terminology. 

The configuration at time t is denoted by a{t) := {ax{t),x G 7?). Let 
Ht denote the distribution of at when we start at time with all vertices 
occupied. 

It is well-known (from a standard coupling argument) that is stochas- 
tically dominated by /ij, if s < t. Hence /xj converges weakly to a probability 
measure, denoted by v (=z^a) as t ^ oo. This measuer v is called the upper 
invariant measure (abbreviated as u.i.m.). Realisations are typically denoted 
by cr = {ax, x G 1?). The occupied cluster of a vertex x (that is, the maximal 
connected component which contains x and of which every vertex is occupied) 
is denoted by C^. (If x is the origin, 0, we often omit the subscript). 

In this paper we study the sizes of occupied clusters under the measure 
v. Motivation comes from work by Liggett and Steif ([25]) who showed that 
for A sufficiently large percolation occurs (that is, z>'a(|C| = oo) > 0), and 
from work by biologists and agricultural researchers. In this latter work (see 
[T9]). limit distributions of contact-like processes (more complicated than 
the 'basic process' described above) were used to model vegetation patterns 
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in arid regions in Spain and North- Africa. In this 'agricultural' context an 
occupied cluster is interpreted as a 'vegetation patch'. For some of these 
models it was claimed in [19j that simulations suggest power law behaviour 
of the cluster size distribution in an interval of some parameter. 

In ordinary percolation models it is known that below the percolation 
threshold the distribution of the cluster size has exponential decay, and that 
power-law behaviour can only occur at the percolation threshold. Triggered 
by the above mentioned claim by biologists and agricultural researchers con- 
cerning very different behaviour in 'their' contact-like processes, we study 
this question for z/^. 

The proof of exponential decay for ordinary (independent) two-dimensional 
percolation goes back to the celebrated paper [20] by Kesten. A crucial step 
in that paper is, somewhat informally and in 'modern' terminology, that if 
the probability of the event A that there is an occupied crossing of a given, 
large, box (square) is neither close to nor close to 1, the expected number 
of so-called pivotal vertices (or, for bond percolation, pivotal edges) is large. 
(These are vertices with the property that flipping the state of the vertex, 
flips the occurrence/non-occurrence of the event A). This step was proved 
in a 'constructive' way, with a 'geometric' flavour. The above mentioned 
large expectation of pivotal vertices implies that the derivative (w.r.t. the 
parameter p) of the probability of A is large. Hence, once the probability of 
A is not very small, a small increase of p makes it close to 1. (This property 
would now be called a 'sharp-threshold' phenomenon'). 

Moreover, by separate arguments, so-called finite-size criteria hold: if the 
probability of A is smaller than some absolute constant e, the cluster size is 
finite a.s. (and its distribution has exponential decay), while if it is larger 
than 1 — e the system percolates. Combining these things gives exponential 
decay of the cluster size for all p smaller than pc- 

Russo (|28]) proved a very general 'approximate zero-one law' and showed 
that the above mentioned sharp-threshold phenomenon can be obtained from 
this more general law, using only a minimum of percolation arguments. In 
particular, in this way Kesten's 'constructive, geometric' arguments could 
be avoided, which is very useful because carrying out such arguments turns 
out to be (too) hard in many dependent models. (We should note, however, 
that for independent percolation the 'constructive' argument still gives the 
shortest self-contained proof, and that in some dependent models, see [2], it 
gives the only currently known proof). 
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Unfortunately, the above-mentioned finite-size criteria involved a so-called 
RSW result of which no ('reasonably general') extension to dependent models 
was known. This explains why, for a long time, Russo's approximate zero- 
one law did not receive much attention in the percolation community. In the 
meantime, results related to Russo's approximate zero-one law, but consider- 
ably sharper and more explicit, were obtained (in other areas of probability 
and mathematics in general) by Kahn, Kalai and Linial ([IB]), Talagrand 
([29]) and Friedgut and Kalai ([Tl]). (See also [27]). 

The importance for percolation of threse sharp-theshold results became 
clear much later, when BoUobas and Riordan ([8]) proved a more robust 
version of the RSW theorem which, combined with a clever use of the sharp- 
theshold results, led to the proof of the long-standing conjecture that the 
critical probability for random Voronoi percolation in the plane is 1/2 (and 
that below 1/2 this model has exponential decay). The robustness of these 
arguments led to similar results for several other two-dimensional percolation 
models (see |9], [11] and [3]). 

The last mentioned paper proved, for 2D lattice models, exponential de- 
cay below the percolation threshold under the quite general condition that, 
informally speaking, the model has a 'nice finitary representation' (in a well- 
defined sense) in terms of finite-valued independent random variables. It 
turned out that under that condition only a weak (not explicitly quantita- 
tive) form, close to that of Russo's ([2H]), of the sharp-threshold results was 
needed. As an example it was shown that the Ising model (with fixed (5 < (5^ 
and external field parameter h playing the role of p in ordinary percolation) 
belongs to this class (thus giving an alternative, more streamlined, proof of 
the main result in Higuchi's paper jlTj). Here the role of finite- valued in- 
dependent random variables was played by the 'independent updates' in a 
suitable discrete-time dynamics. Such a dynamics was possible by (among 
other things) the nearest-neighbour Gibbs property of the Ising model. This 
is a big difference with the contact process, for which we don't know a suitable 
discrete-time dynamics. Therefore, we were not able to derive exponential 
decay for this model from Theorem 2.2 in [3], but instead exploited the full 
quantitative nature of the sharp-threshold results from [T8] and [29j and fol- 
lowed more closely the route used in [8J and [ITj for the Voronoi model and 
the Johnson-Mehl model (which, like the Voronoi model, is a model of planar 
tessellations, but more complicated than te Voronoi model). (Yet another 
route, namely by using results in [15], might work if v would satisfy the strong 
FKG condition, which however (as has been shown by Liggett) it does not. 



4 



We should also note here that the exponential-decay arguments in pQ and 
[26\, which for ordinary percolation work in all dimensions, so far have (even 
in 2D) no suitable analog for dependent percolation). 
Our main result is the following: 

Theorem 1.1. Let A be such that 

z7a(|C| = oo) = 0. 

Then, for every A' < A there exist Ci, C2 > such that for all n > 1 

i^A'dCI > ri) < Ciexp(-C2n). (1) 

Section 2 states properties of the contact process and other (more general) 
ingredients needed in the proof. It also indicates (see the Remark below 
Lemma 12. 3p an alternative proof of the earlier mentioned result by Liggett 
and Steif that percolation occurs for A large enough. 

The proof of Theorem 11.11 is given in Section 3. As mentioned before, 
the essence is still (as it was in ^20|) to show sharp-threshold behaviour for 
certain crossing probabilities. To do this, we follow the main strategy in 
[8J and [llj. However, the model-specific properties of the contact process 
lead to many non-trivial differences in the steps. Therefore, and because the 
contact process is one of the main random spatial models, the proof is given 
in detail. 

We use several well-known results, techniques and terminology from per- 
colation theory. For an introduction to, and general information on, perco- 
lation see [Hj and [ID], and contact processes see 

2 Preliminaries 

2.1 Contact process ingredients 

A well-known classical result for the Contact Process is that there is a critical 
value Ac such that: 

(a) If A < Ac the contact process 'dies out', and u is concentrated on the 
trivial configuration where all vertices are vacant. 

(b) If A > Ac, is non-trivial, and fit converges exponentially to as t ^ 00 
(see Theorem 2.30 and equation (2.31) in [23] which are based on the work 
by Bezuidenhout-Grimmett ([B] and [7]): 
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For all A > Ac there exist C3, C4 > such that for all t > 

^lt{<yo = 1) - i^{<yo = i) < C3 exp(-C4t). (2) 

Since v is dominated by yU^, statement (b) above implies, by standard 
arguments: 

Theorem 2.1. For all A > Ac there exist C3, C4 > such that for all t > 
and all A GG 1? , 

dv{l^t;A,J^X;A) < |A|C3exp(-C4t), (3) 

where dy denotes variational distance, and fit-\ (and i'x-a) are the restriction 
of fit (resp. ux) to A. 

Remark: It is trivial from the definition of Ac that for A below pc no 
percolation of occupied vertices occurs; that is, i^\{\C\ = 00) = 0. As we 
mentioned in the Introduction, Liggett and Steif showed that if A is 

large enough, percolation does occur. It seems to be widely believed (but no 
proof is known yet) that the critical value for having percolation is strictly 
larger than Ac- 

A well-known and very useful way to describe the contact process is by 
means of a space-time diagram (graphical representation): Consider for each 
vertex v Gl? its 'time axis' {v\ x (—00, cx)), and consider five independent 
Poisson point processes on this time axis: one with rate A for each of the 
four directions (left, right, up, down) in the lattice, and one to indicate a 
transition from 1 to 0. The Poisson processes of the different vertices are 
independent of each other. 

The interpretation of a Poisson point on the time axis of v at time t for 
(say) the direction 'right' is that if v is in state 1 at time t, it 'infects' the 
vertex v + (1,0). (That is, if the latter vertex is not occupied, it becomes 
occupied). To visualise this we draw an arrow from (f,t) to [v + (1,0), t). 
For each of the other three directions we act similarly. The interpretation of 
a Poisson point in the fifth process on the time axis of v at time t is that if 
V is occupied (in state 1) at time it becomes immediately vacant (0). In 
the space-time picture this is marked by the symbol * at (f,t) (see e.g. [21], 
Part I, Section 1). 

An active (space-time) path is a path that is allowed to move upward in 
time along the time axes without hitting * points, and to jump from one time 
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axis to another along (and in the direction of) an arrow. For f , w G and 
s < t we denote by {v, s) {w, t) that there is an active path from (t>, s) to 
(w, t). For the contact process starting at time with every vertex occupied, 
a vertex w is occupied at time t > if and only if (in terms of the above 
mentioned space-time diagram) for some vertex v there is an active path from 
(f , 0) to (w, t). In other words, the joint distribution of the random variables 

I{3v e 1} s.t. (t;,0) ^ (w,t)}, w^l? 

is /Xt- Similarly, v is the joint distribution of the random variables 

/{Vt < 3t; G s.t. [v,t) ^ (w,0)}, w G 1? . 

We will often work with the following 'truncated' random variables. First 
some more notation: The distance between two vertices v = (ii,ji) and 
w = (^2,^2) is defined as max(|ii — i2\, \ji — 72!) and denoted by d{v,w). 
The distance d{V, W) between two subsets V and W of Z^ is defined as 
mm{{d{v,w) : v E V,w E W}. For A C Z^, cta denotes the collection of 
random variables (ct„, f G A); straightforward generalizations of this notation 
will also be used. 
Let 

:= /{3(y, t) with d{x, y) = [^/n\ or t = -y/n s.t. {y, t) (x, 0)}, (4) 

and let denote the joint distribution of the random variables 

cri"\xGZ2. 

It is clear from this definition that if A and A' are two finite subsets of 
Z^, and (i(A,A') > 2y/n, then cr^"'* and aj^"'* are independent. It is also clear 
that a is stochastically dominated by cr'^"\ 

From Theorem 12.11 (and simple estimates concerning the 'spatial spread 
of infection in a limited time interval') it follows that 

VA > A, 3C5, Ce > s.t. VA CC I? (5) 
c?v(aA,ai")) < \K\C,eM-C,n^'^). 

Remark: In this paper we often deal with spatial boxes of length of order 
n, and distances of order n to each other; the (somewhat arbitrary) choice 
of -/n in the definition (jlj) is just one of the many possible choices that are 
convenient in such situations. 
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Lemma 2.2. Let Ai, • ■ ■ A^. be 3nxn-rectangles with the property that d{Ai, Aj) > 
2L\/nJ , I < i < j < k. Further let Ai, A2, ■ ■ ■ Ak he events that are completely 
determined by, and increasing in, the a variables on Ai,A2,- ■■ ,Afc respec- 
tively. Then, for every X > X^, 

k 

1=1 



Proof. The first inequality comes from positive association (Harris-FKG in- 
equality) and the last inequality comes from ([5]). The second inequality and 
the equality follow immediately from the definitions. □ 

Let, for a rectangular box R in the lattice, H{R) denote the event that 
there is an occupied horizontal crossing of R. Further, let H{n, m) denote the 
event that there is an occupied horizontal crossing of the box [0,n] x [0,m]. 
For vertical occupied crossings we use a similar notation, with V instead 
of H. From now on, when we write 'crossing' we always mean 'occupied 
crossing'. 

Lemma 2.3. (Finite-size criterion) 

3e > 0\/X> XcBN yN > N 

the following holds: 
(a) If ux{V{3N, N)) < e, the distribution of \C\ has exponential decay. (7) 
(b) Ifux{H{3N,N)) >l-e, then ux{\C\ = 00) > 0. (8) 

Remark: For our purpose (as will become clear later in this paper) we 
do not need N in the above lemma to be uniform in A if A is bounded away 
from Ac. However, although this is not explicitly stated in the literature but 
pointed out to me by Geoffrey Grimmett (private communication), ([2]) and 
related bounds are (by the nature of their proofs in the literature) uniform 
in A (if A is bounded away from Ac). Now, such uniformity would also give 
uniformity of N (in the sence mentioned above). This then, in turn, would 
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(6) 



\i=l 



n^>f (^.) < n ('^a(A) + C^Sn^expi-CV^)) . 



i=l 



i=l 



clearly give an alternative proof of the earlier mentioned result by Liggett 
and Steif that has percolation if A is large enough: Take some A' > Ac- Fix 
N such that for all A > A' statement (b) in Lemma 12.31 holds. It is easy to 
see that (with fixed), if A > A' is large enough, the condition in statement 
(b) holds; hence i^AdCol = oo) > 0. Since this result is already known, and 
not the main subject of this paper, we do not work out the details of such 
alternative proof. (It should also be noted that Liggett and Steif prove more 
than percolation of u: they show (for large A) domination of (high-density) 
product measures). 

Proof, (of Lemma 12.31) . The analog of part (a) was proved for ordinary 
percolation by Kesten in [21], by a block argument. His proof can be, and 
has been in the literature, easily adapted to dependent models with sufficient 
spatial mixing; for an example see Lemma 3.8 in [3J. The mixing property 
described by ([6]) above is 'more than enough' for this purpose: essential is 
that the 'extra term' (here C^Sn"^ exp{—CQ^/n)) in the factors in the r.h.s. of 
(El) goes to as n — >• cxD. 

The analog of (b) was proved for ordinary percolation in by giving a 
suitable (and now well-known) lower bound for the probability of having a 
horizontal crossing of a 4n x 2n box in terms of the probability of the analo- 
gous event for a 2n x n box. If for some n this probability is sufficiently close 
to 1, one can then iterate this procedure and conclude that the probability, 
say rfc, that there is a crossing of a given 2^"'"^n x 2^n box goes very fast to 
1 as /c — > oo. (So fast that ^^^.(l — Tfe) is finite). By Borel-Cantelli it then 
follows that a.s. there is a K, such that for all odd k > K, there is horizontal 
crossing of the rectangle [0, 2''~^^n\ x [0, 2^n] and for all even k > K there is a 
vertical crossing of [0,2^^^] x [0, 2'^"'"^n]. By pasting together these crossings, 
one gets an infinite occupied path; hence, the system percolates. For de- 
pendent percolation models with sufficiently strong mixing properties simple 
modifications of such arguments can be obtained (and have been obtained 
in the literature). Informally speaking, instead of blowing the rectangles up 
by a factor 2, this is then done by a factor 3 to obtain an extra strip in the 
middle of the next rectangle in order to separate the other two strips so that 
the crossing events of these other two strips are almost independent. See, 
for instance, the proof of Theorem 4.8 in [4j for where this has been 

carried out in detail. In practically the same way this can be carried out in 
our current situation (using Lemma 12.21 above in the same way as Lemma 
2.3 was used in the proof of Theorem 5.1 in [1]). □ 
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The following involves what in the Introduction was called a robust ver- 
sion of RSW. 

Proposition 2.4. Let A > Ac- 

If for some p > 0, limsup z/A(-f^(p^, ?^)) > 0, 

n— >oo 

then 

for all p > 0, limsup u\{H{pn, n)) > 0. 

n— >oo 

Proof. A similar result was proved by Bollobas and Riordan [8] for the ran- 
dom Voronoi model (and slightly modified to the above form in [3]). As 
remarked in ^ (see also the first three paragraphs of Section 3.4 in [3]) it 
holds for many percolation models on Z^, namely those that satisfy: (i) A suf- 
ficiently strong mixing property; (ii) A straightforward 'geometric' condition 
about lattice paths (which enables to paste together paths that cross each 
other); (iii) Positive association; and (iv) The condition that v is invariant 
under the symmetries of 7?. 

Lemma [2.21 above is 'more than needed' for (i), and it is easy to see that 
the probability measures v\, A > Ac, also satisfies the other conditions. □ 



2.2 Influence and sharp-threshold results 

Let Vt = {0, 1}" and let Pp denote the product measure with parameter p on 
Q. Let A be an event (i.e. a subset of Q) and let, for 1 < i < n, Ii denote 
the probability that i is pivotal for A. (It is often called the influence of i). 
More precisely, 

Jj := Pp{{uj G Q : exactly one of u and cj*-*^ is in A, 

where cu*-*^ is the configuration obtained from u by flipping the ith. component 
of uj. Talagrand ([22], Corollary 1.2; see also [IB] and [H] for strongly related 
results) proved the following theorem. (Note that 'our' Jj differs a factor 1/p 
from the expression fip{Ai) in Talagrand's paper). 



Theorem 2.5. 



Kp\og{2/p) ^^pmaxjr 
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where is a universal constant. 



Remark: If the event A is increasing (that is, its indicator function is a 
coordinatewise non- decreasing function on Q), the Lh.s. of (Q is, according to 
Russo's formula equal to d/dpPp{A). By this it is easy to see that Theorem 
2.51 implies that if, throughout some interval, say (^1,^2)5 niaxj Jj is 'very 
small', and Pp^(A) is 'not small', then Pp^{A) is 'close to 1'. For such reasons 
Theorem 12.51 and related theorems are often indicated by the name 'sharp- 
threshold' results (besides names like 'influence results'). 

Now suppose there are at least m indices i with the property that Jj = 
maxjlj. There are two possibilities, (a) and (b) below: 
(a) maxj/j < If this holds, then, by Theorem 12. 



A>log(2/p) *^logm'- A>log(2/p) » • ^ ' 



for some universal constant K. 
(b) maxj li > ^^1^. Then, trivially, 



Jj > m max Jj > 



logm 
P 



which is larger than or equal to some universal constant times the r.h.s. of 
( ITOll . Hence (by adjusting the value K if needed) the following holds: 

Corollary 2.6. Let m denote the cardinality of {i : = maxjlj}. Then 

, PJA)(1-PJA))^ 

Remark: The case m = n of this Corollary is, essentially, in ^3] where it 
is derived from the results/methods in p8]; the general case, and its deriva- 
tion from Theorem 12.51 was shown to me by Oliver Riordan (private com- 
munication; see also [I2|)- 

We will use a generalization of Theorem 12.51 and Corollary 12.61 as de- 
scribed below. 

Let fl be as before. Let V G {I, ■ ■ ■ ,n} and let < pi < p2 < I. Let 
P2 denote the product measure on Q under which each component with 
index in is 1 with probability pi and each with index in is 1 with 
probability p2. The generalization of Theorem 12.51 is the following: 
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Theorem 2.7. 



where K' is a universal constant. 

Remark: In [12j (see Theorem 5 in [L2l and the discussion below that 
theorem) it is indicated how to prove Theorem 12.71 by modifications of the 
proofs in Talagrand's paper. An alternative way is to start from the special 
case for p = 1/2 of Theorem 12.51 above, and obtain the full case (and its 
generalization where different coordinates may have a different parameter p) 
from that special case by, informally speaking, representing (approximately) 
the toss of a biased coin, by a combination of tosses of several fair coins. 

From Theorem 12. 71 the following Corollary is obtained in exactly the same 
way as Corollary 12.61 was obtained from Theorem 12.51 



Corollary 2.8. Let m denote the cardinality of {i : Jj = maxj 1^}. Th 



en 



K'p2 log(2/pi) 



3 Proof of Theorem 1.1 



Proof. Let Ai > Ac be such that under the cluster size distribution does 
not have exponential decay. Let A2 > Ai. We will show that z^AadCol = 
00) > 0. This will immediately imply Theorem II. 1[ 

Let Ai be as fixed above, and let e and N = N{Xi) be as in Lemma [2.31 
Let Ln denote a specific in x n rectangle; its precise choice does not matter, 
but for later convenience we choose [n, 5n] x [n, 2n]. By Lemma [2.31 we have 
that 

ux-^{V{3n,n)) > i, for all n> N, 

which by Proposition 12.41 implies limsupj^^^uxiH^Ln)) > 0; so there exists 
an e > and a sequence rii, ^2, ■ ■ ■ such that 

vxAH{Ln,) > £, for all i. (11) 
From now on we consider such fixed sequence. 
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In the introduction to the contact process in the beginning of Section 1 
we assumed that the recovery rate is 1. Of course, the contact process with 
infection rate A and recovery rate 6 is simply a time-rescaled version of the 
contact process with infection rate X/6 and recovery rate 1. In particular, 
these two contact processes have exactly the same upper invariant measure. 
For application of the results in Section 12.21 it is more convenient to work 
with one-parameter Poisson processes for which at each site of the lattice 
the total rate of all the Poisson processes is constant, say 1. Therefore we 
consider the contact process with infection rate g/4 and recovery rate 1 — q, 
where now q G (0, 1) is the parameter. Note that in terms of the space-time 
diagram this means that on each time axis we have a marked Poisson point 
process with density 1, and each point corresponds with a |, | or 

* with probability g/4, g/4, g/4, g/4 and I — q respectively. With respect 
to this new parameter q we use the notation Vq for the law governing the 
above marked Poisson point process, and the notation z/<g> for the upper 
invariant measure of the corresponding contact process. From the above it 
is immediate that 

z><,> = z7_^, gG(0,l), (12) 
or, equivalently, z/a = '^<4A/(i+4A)>, for A G (0, oo). In particular, by (fTTl) . 

z/<,,>(if(L„J) for alH. (13) 
Let jy^g^ be the distribution of {ri^\ x G Z'^) defined by (compare with 

my- 



v!^^ := /{3(y,t) with dix.y) = [v^J or t < -v^ s.t. {y,t) ^"'^'^ (x,0)}, 

(14) 

where {y,t) '•^'.L.^^ {x,0) denotes that there is a space-time path from {y,t) 
to {x, 0) in the space-time diagram with Poisson intensity g/4 for each of the 
four types of arrows, and Poisson intensity 1 — g for *'s. 
It is clear that z/<'^> dominates i^<q>; hence by (|T3l) . 

i?^^^UH{LnJ) > e, foralH. (15) 

Although z/|"^ is, of course, not the same as '^<'4a/(i+4a)>' straightfor- 
ward to get analogs of the earlier 'approximation lemma's'. In particular we 
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get, as analog of ([S]): 



Vg > 4Ae/(l + 4Ae) 3C7, Cg > s.t. VA CC (16) 
dv ('><,>;A,^>2>;a) < |A|C7exp(-C8ni/2). 

A step towards application of the results in Section 12.21 is a suitable 'time- 
discretized' version of z/<'^>. A significant obstacle is to obtain an analog of 
f|T5|) for these discrete variables. 

Recall (see the beginning of this Section) that Ln is the box [n, 5n] x 
[n, 2?T,]. To 'get ample room for the underlying Poisson points' we also con- 
sider the larger box _B„ := [0, 6n] x [0, 3n\. Let z/<"g> be as before. Note that 

the collection of random variables {'ri^\ x E Ln) is completely determined by 
the (marked) Poisson points in the space-time area ST{n) := Bn x [— ra, 0]. 
(In fact only a subset of that area is involved, but for convenience we consider 
this whole area). Let, as before, Vq denote the probability measure govern- 
ing the marked Poisson points. Let a > 0. (Later we choose a sufficiently 
small). Let 6 = n~". We say that an active space-time path vr is 5-stable if 
there are no {y, s) E tt and {z, -u) G x R with d{y, z) < 1 and |m — s| < 5 
such that there is an arrow at {y, s), and a * at {z,u). 

The following Lemma (and the global structure of its proof) is the analog 
of Theorem 6.1 for the Voronoi model in |8j and Theorem 8 for the Johnson- 
Mehl tessellations in [TT] (see also [12] )■ Since the proof is subtle and differs 
in many details from that in [8] and [12] we give a full proof. 

Remark: In some sense the proof of Lemma 13.11 is easier (and shorter) 
than that of the corresponding results in [8] and [11]. This is partly due 
to the fact that in our model the continuous object that has to be properly 
discretized (the time axis) is one- dimensional. This enables us to 'play' with 
the order (in time) of the Poisson points. On the other hand, our model has 
some extra complications, e.g. there is no natural order on the arrow values 
assigned to the Poisson points (an arrow to the right is not always better 
than an arrow to the left). Fortunately these issues can be handled quite 
smoothly. 

Lemma 3.1. (Stability coupling) Let < q < q' < 1. For each n there 
is a coupling ofVq and Vq' such that w.h.p. (i.e. with probability tending to 
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1 as n ^ oo) the following holds: For every x & Ln that has r]x — 1 in the 
first copy, there is a {y,t) G x (— oo,0) with d{x,y) = l^/n\ or t = —y/n 
such that there is a 6-stable space-time path in the second copy from {y,t) to 
{x,0) (and hence 'qx'^ also equals 1 in the second copy). 

Proof. Let ^1 = (So 6 « 5i). We partition every 'time axis' {x} x 

[— oo, 0], x G Z^, in intervals {x} x { — {k + — A; = 0, 1, ■ ■ ■ of length 
5i. From now on, when we use the word 'interval, we will always mean an 
interval of the above form, with x E and {k + l)5i < n. Note that the 
total number of intervals is M„ := |S„|[n/5iJ. Let In denote the union of 
these intervals. 

Note that the total number of Poisson points in X„ is Poisson distributed 
with mean 5iM„. To construct the coupling, first draw a number according 
to the above mentioned Poisson distribution. Now assign N points (called 
'particles') randomly (uniformly, and independently of each other) to the 
above mentioned set X„. If a particle is assigned to the space-time location 
{x, t), we say that its time coordinate is t. Call an interval 'occupied' if it has 
at least one particle. Call two different intervals {x} x {—{k+ l)6i, —k6i] and 
{y} X (— (/ + l)5i, —ISi] neighbours if d{x,y) < 1 and \k — l\ < 1. This gives 
rise, in an obvious way, to the notion of 'clusters (of occupied intervals)'. 
(This notion of cluster is of course different from that introduced earlier in 
this paper. Since this 'new' notion of cluster is used only in this proof (and 
the other notion is not used here), this should not cause any confusion and 
we even use the same notation C). 

We have already assigned, to each particle, a precise location in X^- How- 
ever, we 'suppress' this precise information and only 'keep' the following 
partial information: For each interval the number of particles assigned to it, 
and for each occupied cluster of intervals the relative order (w.r.t. their time 
coordinates) of all particles in that cluster. 

Prom now on, when we mention a cluster C, we mean not only its corre- 
sponding set of intervals, but also the above mentioned partial information. 

By the size of a cluster we mean the number of particles in the cluster. 
By a simple 'counting argument' it is easy to see that there is 7 > (which 
is a function of 5i which goes to 00 as (5i — >^ 0) such that the probabihty that 
the size of the cluster containing a given interval is larger than k is at most 
exp(— 7A;), and hence 
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P(3 a cluster with size > Clogn) < M„(5i) exp(-7Clogn) < C'n^+^^Z^-^^, 

(17) 

where C is a constant. Since, for fixed a, 5\ is a function of n, we can 
consider 7 as a function of n satisfying 7(n) ^ 00 as n ^ 00. Hence, for 
each choice of a, we have that for every /3 > 0, w.h.p. (i.e. with probabihty 
tending to 1 as n — > 00) all clusters are smaller than /3 log n. 

Let C be a cluster (including the partial information given above). Now 
consider, for both parameter values {q and gr') the conditional distribution 
of the precise configuration for C, i.e. the types (*, , | or [) and 
precise locations of all particles in C, given the partial information. The two 
conditional distributions can be coupled by the following natural procedure 
which gives two 'typical realizations' (one for each of the two parameter 
values) of the precise configuration. 

The first step in this procedure is to assign, with equal probabilities (1/4), 
a tentative f or | to each particle (independent of the other particles) 

of C. The interpretation is that eventually a particle is chosen to represent 
an arrow, the type of arrow is exactly its above mentioned tentative one. We 
also assign to each particle i, independent of the other particles, a random 
variable Ui, uniformly distributed on (0,1). These variables will be used 
below to decide if a particle corresponds with an arrow or with a *. 

The next step is to go from relative order of positions to precise positions. 
Consider the conditional distribution of the precise time coordinates of the 
particles of C, given their (already known) relative order in time and the 
intervals they are located in. Now simply assign the precise locations by 
drawing from this distribution. Later we will refer to this procedure as the 
'time assignment procedure'. 

Note that both steps above are the same for both 'reahzations' (the one 
for parameter q and the one with parameter q'). However, the next (and 
final) step, in which the types of the particles are fully determined, will take 
into account the parameter value: For each particle i of C do the following. 
If Ui < q, the type of i (in each of the two copies) is equal to the earlier 
mentioned tentative arrow. If Ui e {q,q'), its type is * in copy 1 and equal 
to the tentative arrow in copy 2. If Ui > q', the type is * in both copies. 

Now we have two realizations, say ujc{l) and a;c(2), and it is easy to see 
that they are 'typical' w.r.t. the two conditional distributions mentioned 
above, so that we indeed have a coupling of these two conditional distribu- 
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tions. Also note that co'c(2) > uJc{l) in the sense that each particle in co'c(l) 
that has an arrow-type, has the same arrow-type (and space-time location) 
in co'c(2). Let this coupling be denoted by Pc- 

Doing this for each cluster, independently of the other clusters, gives a 
natural coupling of the two probability measures in the statement of the 
lemma. However, it is not yet what we want: Although it satisfies the 
property between brackets at the end of the Lemma, it does not necessarily 
satisfy the stability property in the lemma. The coupling we do want is 
obtained as follows, where we go back to the level of a given cluster C. Recall 
the two copies ujc{^) and ujc{2) above, and their joint distribution Pc- From 
Pc we will construct a modified distribution Pc of which the two marginal 
distributions are the same as those of Pc- To avoid an abundance of notation, 
we will drop the subscript C from uc{l) and ujc{2). 

Recall the time assignment procedure in the (second step of the construc- 
tion of) Pc- Let B be the event that in uj{l) (and hence, since the particle 
locations for uj{l) and uj{2) are the same, also in uj{2)) there are two differ- 
ent particles in C whose time coordinates differ less than 6. The probability 
of B (or, more precisely, the conditional probability of B given the partial 
information on C) is maximal if C consists of one interval only, in which case 
it is less than or equal to \Cf^ = |C|^2n~°/^, where \C\ denotes the number 
of particles in C; so 

Pc{B) < |C|^2n-°/^ (18) 

Recall the use of the variables Ui in the determination of the types of the 
points. Let G be the event that each particle in a;(l) is of type * and each 
particle in a; (2) has an arrow type. Note that this event happens if and only 
if Ui e {q, q') for all particles i in C, so that we have 

Pc{G) = {q'~qf\. 
By this and (ITS]) there is a /?' > such that 

Pc{G)>Pc{B)ii\C\<P'\ogn, (19) 

and n is sufficiently large. 

Now let B' denote B\G. If \C\ < /3'\ogn then, by 1^ there is a 
(measurable) subset G' C G \ B and a 1-1 map ip : B' ^ G' with the 
property that ip and ip'^ are Pc-preserving. To each pair {ijj{1),uj{2)) G B' 
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this map assigns the pair 

V^(^(l),a.(2)) = ((V^(^(l),.;(2)))(l),(V^(^(l),.;(2)))(2)). 

Now a modified couphng, called Pc, is obtained from Pc by exchange, between 
B' and G', of the second copy, using ip, as follows. (Such a modification 
is called a 'cross-over' in [8]). If \C\ > f3'\ogn we simply take Pc = Pc- 
Otherwise, a typical pair (a;(l),a)(2) under Pc is drawn as follows. First draw 
a pair tu(2)) under Pc. If (tu(l), cj(2)) G {B' U G")'^ take (tu(l), tD(2)) 

equal to (cj(1), cj(2)). If cu(2)) G fi', take (::;(!) = to and (:D(2) = 

(?/^(cu(l),cj(2)))(2). Finally, if cu(2)) G G" take cu(l) = cu(l) and (1^(2) = 

co'(2)))(2). Since, in all cases, (a'(l) = uj{2), it is immediate that 
the first marginal of Pc is equal to that of Pc- A short inspection shows that 
also the second marginal of Pc is equal to that of Pc- 

Now the 'overall' coupling of Vg and Vg' announced in the statement of 
the Lemma is obtained, in a natural and straighforward way, by constructing 
the pair (a)c(l), u;c(2) for each cluster C separately, independently of the other 
clusters. 

To check the required properties of this coupling, first look again at one 
single cluster C. Suppose that \C\ < P'logn. Let {uj{1),lj{2)) and the corre- 
sponding pair (a)(1), u;(2)) be as above (so, in particular, u;(l) = uj{l)), and 
suppose that a)(l) has a certain space-time path vr, with at least one arrow, 
within C. We have to consider a few cases: 

(i) If {uj{l),uj{2)) G {B' U G'Y then cD(2) = uj{2) > uj{l) = Cu{l). Hence vr is 
also an active space-time path for u'. Now distinguish two subcases: 

(i(a)) If (ti;(l), a;(2)) ^ B then, by the definition of B (and noting that the 
time coordinates of a particle in C and a particle outside C differ at least 
6i > 6) 71 is 5-stable. 

(i(b)) If (w(l),w(2)) G G, then (by definition) cu(2)(= lu{2)) has no *'s in 
C; hence (see the note in (ia) vr is 5-stable. 

(ii) If (ci;(l), a;(2)) G B', then ip{u{l),u{2)) G G, so each particle in the con- 
figuration a)(2) represents an arrow. In particular u)(2) has all the particles, 
and their corresponding arrows, of vr. Although the precise locations of these 
particles may have changed a bit, their relative order is the same, and is easy 
to see (see again the note in subcase (ia) and use that there are no * particles 
in Cj{2)) that a;(2) has a 5-stable path within C from the the beginning of tt 
to the end of tt. 

(iii) Finally, observe that {uj{1),uj{2)) can not be in G', because of the men- 
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tioned space-time path vr in u{l) (= uj{l)). 

Concluding, the only way in which oj{2) has no 5-stable path from the 
the beginning of vr to the end of vr is when \C\ is larger than /?' logn. 

Using the above mentioned property of the single-cluster couplings (and 
noting that the condition, in Lemma I3.H that r^i = 1 in the first copy, im- 
plies the existence of an active spacetime path to {x, 0) in the first copy from 
some {y,t) satisfying d{x,y) = [y/n\ or t = —^/n), yields a similar property 
for the 'overall' coupling of Vq and Vq'. This, together with the fact that, 
as we saw before (see f|T7|) ). the probability that there is a cluster C with 
\C\ > P'logn goes to as n — ^ oo, completes the proof of Lemma \3A] 

□ 

Fix a value q in the interval {qi,q2), where qi = 4Ai/(4Ai + 1) as before 
(see (US])) and = 4A2/(4A2 + 1). 

Now we are ready to introduce — 1 valued random variables to which 
we can apply the results in Section I2.2[ Let the box Bn and the space-time 
region ST{n) be as before (see the paragraph preceding Lemma 3.1). Now 
partition every time axis in intervals of length 6, with S as defined a few lines 
before Lemma [3.11 

As before, we have on each time Poisson point process with density 

1, and each Poisson point is, independently of the others, of type * with 
probability 1 — g, and of each of the types — , t, | with probability g/4. 
Let f G -Bn and k E N, < k < n/6. By the kth interval of v (for the above 
mentioned partition) we will mean {v} x (— M/4, (— A; + 1)5/4], and we define 

j^iv,k,s) ._ j|g ^ Poisson point of type * in the kth interval of f }. 
Similarly define 

poin of type in the kth interval of v}, 

and, analogously, X^'^'^\ X^^'^'^^ and xj^''^'^''. Note that this is a collec- 
tion of independent — 1 valued random variables. 

Recall the definition of rji^^ and jy^q^ below equation f|T3l) . The X vari- 
ables defined above, give only 'crude' information about the space-time di- 
agram: they tell which of the types *, — > etc. occur in each interval, but 
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they do not tell their precise locations inside the intervals. Nevertheless, this 
incomplete information is often enough to conclude that there is a certain 
space-time path. Let t]^'^^ be the indicator of the event that the values of 
the X^''''^^ variables imply that r/l,"'' = 1. 

Remark: Note that if r]^'^ = 1 then, after for some Poisson points with 
mark * this mark is replaced by an arrow, still r]i"-^ = 1. The same remark 
holds for rji^'^^ instead of rii"'\ 
It is easy to see that 



vi""^ > vi""'^^ > t) with d{v, w) = [Vn\ or t = - s.t. (20) 

3 5-stable space-time path from {w,t) to (f,0)}. 

Hence, with the following notation (where i? is a box) 

H^"''^\R) := {3 //^"'■'^-occupied horizontal crossing of R}, 

we get 



Vq (i/("*'^)(L„J) > Pg,(3r/("»)-occupied hor. crossing of L„J - o{ni) 
= i.(r;^>(/7(L„J-o(n,), 

where the inequality comes from the second inequality in (120|) and Lemma 
13.11 and the equality comes directly from the definitions. 

By f fTSj) . and obvious monotonicity (see the Remark preceding fl20|) ). this 
gives the following Lemma: 

Lemma 3.2. For each choice of a the following holds for all sufficiently large 
i: 

P,((i/(-'^)(L„J) > |, q>q. (21) 

Now we "wrap around the box horizontally" (by identifying every 
vertex {6ni,y) on i?„. with the vertex {0,y)), thus turning this box into a 
cylinder. Define rji""^''"'' as the natural analog, for the cylinder, of r]i^"^\ 
Remark: 'By the truncation to distance ^/n of these variables', and because 
the left and right side of L„- have distance larger than ^/n to the boundary 
of Bn^, the event that there is a r^.*-"'"'^''" ''-occupied horizontal crossing of L„. 



20 



and the event that there is a rj -occupied horizontal crossing of are 
the same. 

Let A*^""^) be the event that at least one of the (Grij — 1) horizontal translates 
of L„. on this cylinder has an ?7f"'*'^''^''-occupied horizontal crossing. Note that 
j\{ni,&) jg g^jjj defined in terms of the random variables X ''^'^ defined earlier. 
For each choice of a the following holds for all sufficiently large i: 

> (i7K.^)(L„J) > |, g > g, (22) 

where the first inequality is (taking into account the above Remark) trivial 
and the last inequality is exactly Lemma 13. 2[ 

As said before, the X variables are independent — 1 valued random 
variables. Further, for each v and fc, x'f'^'^'^ has probability 1— exp(— (1 — g)(5) 
to be 1. If z is large enough (so that 5 is small) this is larger than or equal 
to (1 — q)6/2. Each random variable I^'^'^') has probability 1 — exp(— (5g/4) 
to be 1, which is larger than q6/8 for i sufficiently large. The same holds for 
the other three arrow types. 

Also note that the event A^"^'^^ is partially symmetric in the following 
sense: For fixed value k and fixed < / < 3n all variables X^^'^''^ with v G Bn 
with y— coordinate /, 'play the same role'. In particular, each of them has 
the same probability to be pivotal for the event A^'^^'^\ The same statement 
holds for each of the other three arrow types and for type *. Further note 
that for each k and / the number of such random variables X^''''^^ is of order 
n. (And, again, the same statement holds for each of the other types). By 
Corollary 12.81 (with m equal to our 'current' n) and a simple modification of 
Russo's formula this yields 



^^^^^ ^ M275^ ' ^ ^ ^^^^ 

where Cg > depends on q and q2 only. 

Let £* > 0. By (1221) . ( l23l) and because Vq{A^'^^'^'>) is clearly non- decreasing 
in g, it follows that, for every choice of a, the following holds for all sufficiently 
large i: 

If P,2(A("-^)) < 1 - £*, then, for all g G [g, ga), 

dq^'^^^ ^-^^2^ log(2/5) - a ' 
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(where the last inequahty used that 5 = n and hence 

By choosing a sufficiently small this gives the following Lemma: 

Lemma 3.3. For every e* > there is an a > such that for all sufficiently 
large i 

Pg^l^^""''^) > 1-^*- (24) 

Now, if there is a horizontal crossing of one of the above mentioned trans- 
lates of Lni, there must be a horizontal crossing in the 'hard' direction of at 
least one of the following (six) translates (on the cylinder) of the rectangle 
[0,3ni] X [ni,2ni]: 

[jni, {j + 3)ni{ mod Gui)] x [m, 2ni], < j < 5. 
Hence, by the usual 'square root trick', 

P,,(i7("-^)([0,3n,] X K,2n,]) > 1 - {l-V,M^--'^f\ 

which, combined with Lemma [373| immediately gives that for every e* > 
there is an a > s.t. for all sufficiently large i 

V,, (/7("-^)([0,3n,] X [0,n,])) >l-6*. (25) 
Finally, the following Proposition is obtained: 
Proposition 3.4. 

lim ux2iH{3ni,ni) = 1. (26) 

Proof. Let > be given. By (fT2!) (and the definition of ^2) i^\2 {H{3ni, Ui)) = 
i/<g2>(-ff(3?2j, nj)). Hence, by f|T6l) . liminfj^oo i^A2(-f^(3ni, n^)) is equal to 

lim infj^oo i^<^2>(-^(3'^i; ■"-«)) which, by the ffist inequality in (12U]) . is larger 
than or equal to 

liminfP,, (i7("-'^)([0,3ni] x [0,ni])) . 

This last expression is, by (the statement ending with) fl25|) and a suitable 
choice of a, larger than 1 — e* . Summarizing, we have that for every e* > 0, 
liminfj^oo ^\2{H{3ni,ni)) is larger than 1 — e*. □ 
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Proposition 13.41 together with the finite-size criterion Lemma 12.31 imme- 
diately yields z>'a2(|Co| = oo) > 0, which, as observed in the beginning of this 
Section, completes the proof of Theorem 11.11 

□ 
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